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Ocular dominance patterns in mammalian visual cortex: A wire length 
minimization approach 

Dmitri B. Chklovskii and Alexei A. Koulakov 

Sloan Center for Theoretical Neurobiology, The Salk Institute, La Jolla, CA 92037 

We propose a theory for ocular dominance (OD) patterns in mammalian primary visual cortex. 
This theory is based on the premise that OD pattern is an adaptation to minimize the length of 
intra-cortical wiring. Thus we can understand the existing OD patterns by solving a wire length 
minimization problem. We divide all the neurons into two classes: left-eye dominated and right-eye 
dominated. We find that segregation of neurons into monocular regions reduces wire length if the 
number of connections with the neurons of the same class differs from that with the other class. 
The shape of the regions depends on the relative fraction of neurons in the two classes. If the 
numbers are close we find that the optimal OD pattern consists of interdigitating stripes. If one 
class is less numerous than the other, the optimal OD pattern consists of patches of the first class 
neurons in the sea of the other class neurons. We predict the transition from stripes to patches 
when the fraction of neurons dominated by the ipsilateral eye is about 40%. This prediction agrees 
with the data in macaque and Cebus monkeys. This theory can be applied to other binary cortical 
systems. 



I. INTRODUCTION 

In the primary visual area (VI) of many mammals, 
most neurons respond to the stimulation of the two eyes 
unevenly: they are either left-eye or right-eye dominated. 
In some species, left- /right- eye dominated neurons are 
uniformly intermixed in space. In others, left-/right-eye 
dominated neurons are segregated resulting in a system of 
alternating monocular regions. This system is known as 



the o cular dominance (OD) pattern ( Wiesel and Hubcl, 
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Moat thcoriota interested in the OD pattern ( Erwin ct 
al., 1995), ( Swindale, 1996| ) have been modeling its devel- 
opment. They succeeded in generating OD patterns of re- 
alistic appearance. However, several why rather than how 
questions remained unanswered. Why, from the func- 
tional point of view, do the OD patterns exist? Why do 
some mammalian species have OD patterns while others 
do not? Why do monocular regions have different ap- 
pearances (stripes as opposed to patches) between differ- 
ent species and even between different parts of VI within 
the same animal? 

Mitchison, 1991, suggested an answer to the first ques- 
tion using the wiring economy principle (|Cajal, 1995), 



(Airman and Kaas, 1974), (Cowey, 19791), (Chcrniak 
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1992|), flYoung, 1992| ), flChklovskii and Stevens, 1999| ) 
The idea is that the evolutionary pressure to keep the 
brain volume to a minimum requires making the wiring 
(axons and dendrites) as short as possible, while main- 
taining function. In general, the function of a cortical 
circuit specifies the connections between neurons (wiring 
rules). Therefore the problem presented by the wiring 
economy principle is to find, for given wiring rules, the 
spatial layout of neurons that minimizes wire length. 
Then we can understand the existing layout of neurons 
as a solution to the wire length minimization problem. 



We adopt the wiring economy principle and address the 
above questions by formulating and solving a wire length 
minimization problem. Because of the columnar organi- 
zation of the cortex (Mountcastle, 1957) we consider a 
two-dimensional neuronal layer of uniform density. The 
number of the left-eye dominated neurons is a fraction 
jjj of the total number, and fn is a fraction of right-eye 
dominated neurons (//, + fn. = I). 

We consider only intra-cortical connections be- 
cause they constitute the majority of gray matter 
wirin g (|LeVay and Gilbert, 1976|), (Peters and Payne. 



1993| ) 



t eVay and Gilbert 197j 
hmed et al., 1994 ) allowing us to neglect the 



thalamic afferents and other extra-cortical projections. 
We assume that each neuron receives synapses from N s 
neurons dominated by the same eye and from N neurons 
dominated by the opposite eye. In other words, because 
synapses are unidirectional the resulting wiring rules re- 
quire each neuron to get unidirectional connections from 
N s neurons dominated by the same eye and from N neu- 
rons dominated by the opposite eye. 

Given these wiring rules we look for an optimal layout 
of neurons which minimizes the total length of connec- 
tions. Depending on the values of N Sl N , and /l, opti- 
mal layout belongs to the one of the four phases shown 
in Fig.|lj where left-eye dominated neurons are shown in 
black and right-eye dominated neurons - in white. In the 
Salt and Pepper phase left-eye and right-eye dominated 
neurons are uniformly intermixed, Fig.|lja. The Stripe 
phase consists of alternating monocular stripes of neu- 
rons dominated by either eye, Figjljb. The L-Patch con- 
sists of the patches of the left-eye dominated neurons sur- 
rounded by the right-eye dominated neurons, Fig.plc. The 
R-Patch consists of the patches of the right-eye domi- 
nated neurons surrounded by the left-eye dominated neu- 
rons, FigJ^d. 

Our approach differs from Mitchison's in that we drop 



the retinotopy requirement, that is our wiring rules do 
not take into account receptive field positions. This sim- 
plification is supported by the existence of the receptive 



field scatter (Hubel and Wiesel, 1974), random variation 
in the receptive field position between adjacent neurons. 
For example, in a macaque retinotopy exists only on the 
scales greater than « 1mm (Hubcl and Wiesel, 1974), 
which exceeds the typical size of monocular regions. By 
disentangling the retinotopy from the OD problem we 
simplify it enough to map out a phase diagram. 

In the Discussion we compare our predictions with 
the data from macaque and Cebus monkeys and find 
good agreement. Also, we discuss simplifying assump- 
tions made in the paper and possible ways to extend the 
theory. 
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FIG. 1. Different appearances of the ocular dominance pat- 
tern. Left-eye dominated neurons are black while right-eye 
dominated neurons are white, (a) Salt and Pepper phase, uni- 
formly intermixed left/right neurons. (b) Stripe phase, alter- 
nating monocular stripes, (c) L-Patch phase, circular left-eye 
islands in the right-eye sea. (d) R-Patch phase, circular right- 
eye islands in the left-eye sea. 



dictions of the critical value agree with the data from 
macaque and Cebus monkeys. 




FIG. 2. Ocular dominance phase diagram calculated in the 
lattice model. Optimal phases are shown as a function of 
the relative number of the same-eye to other-eye connections 
N s /N and a fraction of left-eye neurons /. Range of the 
Stripe phase optimality is shown in black, Patch phase - in 
white, Salt and Pepper - in grey. 

In the following section we formulate a wiring problem 
on a lattice. For small N s and N we solve it analytically 
while for large N s and N we solve it numerically. Results 
are shown in Fig.|2|. Next, we introduce a continuous for- 
mulation of the problem. We prove that Salt and Pepper 
is an optimal layout when N s = N . Then we show that 
for N s ^ N segregation of neurons according to their 
OD reduces wire length. We calculate in perturbation 
theory the wire length for Stripe and Patch phases and 
find the range of parameters for the optimality of each 
phase. Perturbation theory provides an analytical treat- 
ment of neuronal clustering so common throughout the 
nervous system. The calculated phase diagram is similar 
to that obtained in the lattice model. 



II. RESULTS 



A. Lattice model 



We present the central results of the paper on a phase 
diagram, Fig.^, showing optimal phases for various ratios 
of same-eye to other-eye connections N s /N and fractions 
of left-eye neurons fz- If the numbers of same-eye and 
other-eye connections are equal, N s /N — 1 then Salt 
and pepper phase is optimal. Otherwise, if N s /N ^= 1 
the wirelength is minimized by an OD pattern consist- 
ing of alternating monocular regions. The shape of these 
regions depends on the relative fraction of the left-eye 
dominated neurons, fz- When the numbers of neurons 
dominated by each eye are close, fz « Jr, the Stripe 
phase is optimal. When the fraction of left-eye (right-eye) 
dominated neurons drops below a critical value f c ~ 0.4 
the L-Patch (R-Patch) phase becomes optimal. Our pre- 



Although the arrangement of neurons in cerebral cor- 
tex is anything but grid-like we can understand many fea- 
tures of the neuronal layout by studying lattice models. 
These models compensate in clarity and computability 
what they lack in realism. Of course, we need to make 
sure that the results are independent of the particular 
choices of lattice parameters (for example the number of 
nearest neighbors). 

We consider arranging a large number of neurons on 
a two-dimensional square lattice. Each site must be oc- 
cupied by either left-eye or right-eye dominated neuron. 
The number of the left-eye dominated neurons is a frac- 
tion fz of the total number of neurons and fn is a frac- 
tion of right-eye dominated neurons. The problem is to 



find a layout which minimizes the total length of wiring 
specified by the following rule. Each left-eye neuron has 
unidirectional connections with N s left-eye neurons and 
with N right-eye neurons. Each right-eye neuron has 
unidirectional connections with N s right-eye neurons and 
with N left-eye neurons. Unidircctionality of connec- 
tions means that connecting neuron A to neuron B, does 
not necessarily imply that neuron B connects to neuron 
A. The motivation for this rule comes from the unidirec- 
tional properties of synapses in the brain. 

Because we attempt to minimize wire length we as- 
sume that for a given layout the connections are estab- 
lished optimally. Thus the problem is reduced to com- 
paring optimal wiring for various layouts. Therefore, we 
will assume that each neuron makes the shortest possible 
connections satisfying wiring rules. 

Small numbers of connections per neuron. We start 
by finding optimal layouts for three illustrative examples 
of wiring rules with small numbers of connections, N s 
and N . We caution the reader that because of the small 
numbers of connections phase assignments may seem ar- 
bitrary. These examples are chosen to illustrate our main 
results which will be confirmed both in the lattice model 
with large N s and iV later in this section and in the 
continuous model (section IIB). 

For the first two examples we set equal numbers of 
left and right neurons, //, = fa = 1/2. In the first ex- 
ample each neuron connects with equal numbers of the 
same-eye and other-eye neurons, N s = N — 2. Then 
the optimal layout is the system of alternating rows of 
left /right neurons, Fig.pk. This layout is a realization of 
the Salt and Pepper phase, Fig.pk, because each neuron 
has an equal number of left and right neurons among its 
nearest neighbors. To calculate the length of connections 
per neuron, I, we notice that in this layout all neurons 
have the same pattern of connections. By considering 
one of them, Fig.pk, we find that I = 4. This layout is 
optimal because each neuron makes all of its connections 
with nearest neighbors. 



Fig.|b, we find I = 3 + V2 « 4.41, greater than I = 4 
for the Salt and Pepper phase. Here each neuron has 
among its nearest neighbors only one other-eye neuron, 
while the wiring rules require connecting with two other- 
eye neurons. A connection to the next nearest neighbor 
is longer making the layout suboptimal. We confirm the 
optimality of the Salt and Pepper phase for N s = N 
both numerically for large N s , N and analytically. 

In the second example each neuron connects with more 
same-eye than other-eye neurons: N s = 3, N = 1. Then 
a realization of the Salt and Pepper phase, Figfik is not 
optimal anymore. The length of connections per neuron 
is I ~ 4.41, while the Stripe phase, Fig.|]b gives I = 4. The 
Salt and Pepper phase loses in wiring efficiency because 
there are not enough same-eye neurons among nearest 
neighbors and connections with the next nearest neigh- 
bors are needed. The Stripe phase, Figjjb rectifies this 
inefficiency by having each neuron make connections only 
with nearest neighbors. Thus, clustering of same-eye neu- 
rons is advantageous if each neuron connects more with 
the same-eye than with the other-eye neurons. 
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FIG. 4. Ocular dominance patterns for fi, = 1/2 and 
N s = 3, N — 1. (a) A realization of the Salt and Pepper 
is suboptimal. (b) A realization of the Stripe phase gives 
minimal wire length. 
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FIG. 3. Ocular dominance patterns for /l = 1/2 and 
N s = N — 2. (a) A realization of the Salt and Pepper phase 
gives minimal wire length, (b) A realization of the Stripe 
phase is suboptimal. 

A suboptimal layout for the same wiring rules is illus- 
trated by a realization of the Stripe phase, Fig.pp. In this 
layout each neuron has the same pattern of connections 
up to a mirror reflection. By considering one of them, 
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FIG. 5. Ocular dominance patterns for fn — 3/4 and 
N 3 — 3, N — 1. Realizations of the Salt and Pepper (a) the 
Stripe (b) are suboptimal. (c) A realization of the L-Patch 
phase gives minimal wire length. 

In the third example we use the same wiring rules 
(N a — 3, N — 1) but take different numbers of left/right 
neurons, //, = 1/4, /# = 3/4. The realizations of the 
Salt and Pepper phase is shown in Fig.pk and of the 



Stripe phase in Fig.gb. In these layouts, different neu- 
rons have different patterns of connections. To find the 
wiring length per neuron we average over different pat- 
terns and find for the Salt and Pepper phase / « 5.02 and 
for the Stripe phase I — 4.5. A more efficient layout is the 
L-Patch phase, Fig.pp, where I ~ 4.21. Although we can- 
not prove that the L-Patch phase is optimal, this seems 
likely. Thus, the optimal shape of monocular regions de- 
pends on the relative numbers of left/right neurons. 

Large numbers of connections per neuron. Lattice 
models with small numbers of connections per neuron 
yield quick results good for illustration purposes. How- 
ever, they are difficult to generalize to the wiring rules 
with large numbers of connections more appropriate 
for cortical circuits where each neuron connects with 
s=s 10 neurons. For example, attributing the layouts in 
Fig.||a,[|a to the Salt and Pepper phase rather than the 
Stripe phase may seem arbitrary. Therefore we study 
lattice problems with large numbers of connections per 
neuron. 

When the number of connections per neuron is much 
greater than the number of nearest neighbors, the effect 
of the discreteness of the lattice on the results is negligi- 
ble. In particular, for a given fraction fi, and ratio N s /N 
only the periodicity of the optimal layout depends on the 
N . Thus the solution of the wire length minimization 
problem for one value of N can be generalized to other 
problems with the same /z, and N s /N . 

We solve the wire length minimization problem using 
the following numerical algorithm. We fix the values of 
Jl and N s /N . We consider neuronal layouts belonging 
to several phases: Salt and Pepper, Stripe, L-(R-) Patch 
(both triangular and square lattice), Checkerboard (only 
f° r /l = 1/2). For each phase we find the optimal pe- 
riod which minimizes wirelcgth. Then we compare wire 
length in the optimal layouts of different phases. We plot 
the optimal phases for various values of Jl and N s /N on 
the phase diagram, Fig. pi These results were discussed 
above. 



B. Continuous model 

In this section we study the limit when N s and N are 
very large. Instead of considering each receiving connec- 
tions neuron separately it makes sense to treat them as 
a mixture of two " liquids" , the left- and right- eye ones, 
having continuous in space densities. Segregation of such 
a mixture implies that the OD structure is formed. 

The model. In this subsection we will assume that the 
neuron configuration represented by the neuron densities 
is given. It can be any arbitrary configuration, including 
Salt and Pepper, Stripes, or Patches. For the given neu- 
ron densities we draw the connections between the cells 
which: 



• Satisfy the wiring rules (N s same and N other neu- 
ron connection have to be established); 

• Minimize the total wire length. 

In the end of the subsection we calculate the wire length 
for the Salt and Pepper configuration. 

We consider the mixture of neurons of two types: domi- 
nated by the left and right eyes. We assume that the neu- 
rons are located in the plane. This assumption is based 
on the fact that the OD remains constant in the direc- 
tion perpendicular to the cortex surface. The variables 
of the problem can therefore be considered functions of 
the remaining two coordinates, r. 

Instead of considering each individual cell we charac- 
terize the neuron configuration by continuous local vari- 
ables. We define the local density of neurons dominated 
by the right eye n R (r) as the average density in a square 
containing sufficiently large number of cells (> 10), yet 
small compared to the typical spatial scales of the con- 
figuration (~ 1mm). Similarly ul (r) is the local density 
of cells dominant by the left eye. Although both n R (r) 
and Ul (r) can vary in space, the total density of neu- 
rons uq = n R (r) + nj, (r) is a constant, independent of 
the position in the cortex. 

In our model n R (r) and nr, (r) completely define the 
neuron configuration. For example the Salt and Pepper 
configuration, in which the densities of right-eye and left- 
eye neurons are uniform, can be defined as follows 



n R (r) = n R = f R n , 
n L (r) eS l = f L n , 



n R + n L — n , 



(1) 



Here f R is defined as the fraction of the right-eye neurons 
with respect to the total number of cells (in general not 

1/2). ^ 

Having defined the neuron configuration by fixing the 
densities n R (r) and ni(r) we proceed to establishing the 
connections between cells. Two requirements have bo te 
taken into account. First, we have to satisfy the wiring 
rules. Second, for given densities n R (r) and iiL(r) the to- 
tal length of connections has to be minimum. Consider 
a pattern of connections from a neuron dominated, for 
example, by the right eye. Consider also the region in 
the cortex it is connected to. There are in fact two such 
regions, for the right- and left-eye connections. We claim 
that each of these regions is a disc. To prove this, notice 
that if the connections are produced with neurons out- 
side of this disk rather than inside the wire length is in- 
creased. This is inconsistent with the requirement of the 
optimum wiring for a given configuration. We denote the 
radii of these two disks R RR (r) and R R l (r) , implying 
the radii of right-eye neuron at point r connection regions 



to the right-eye and left-eye cells correspondingly. Simi- 
lar quantities can be introduced for the left-eye neurons 
at point r, i. e. Rlr (r) and Rll (r). We introduce the 
index notation i — {R, L}. Then the four radii discussed 
can be collapsed into one notation R ik (r), standing for 
the radius of the connection region for the neuron of OD 
i at point r to the cells of OD k. The radius can be de- 
termined from the wiring rules (N s and N connections 
to the cells of the same and other OD respectively have 
to be established): 



N ik = I dr'n k (!•') 

|r-r'|<i^(r) 



(2) 



Here the elements of matrix Nik, i = {R, L} are equal to 
N s if i = k and N otherwise. 

It is now possible to determine the total connection 
length in the cortex C. To this end we add up the lengths 
of the connections of individual neurons Ln~ (r) over the 
whole area: 



•* : i n> T 



i,k=R,L 



where 



Lik (r) 



v-r'\<R lk (r) 



dr'n k (r') |r — r' 



(3) 



(4) 



The last factor in this expression is the connection length 
as a function of separation r — r' between neurons. In 
principle, cost function may not be a linear function of 
separation. However, we take it to be linear for the sake 
of simplicity. Eqs. (0) - (Q) define our model completely. 
Using Eq. (0J we calculate the wire length for the ho- 
mogenous Salt and Pepper configuration. To this end we 
substitute the densities given by Eq. (Il]) into (M) to find: 



SP 



n ik 



TTHk 



Then using Eqs. (0) and M) we obtain 



SP 



Lik - ~R'ik 



(5) 



(6) 



and finally 



C 



SP 



2A 




(7) 



»R 



\fnk~ 



where A is the total area of the cortex. 

In the next subsection we show that wire length can 
be reduced with respect to (R) by introducing a small 



inhomogeneity into the neuron densities ur and ul- To 
this end we treat our model (J^) - (|4J) in the framework of 
the perturbation theory. 

Instability of the Uniform State Leads to the Forma- 
tion of Patterns. The purpose of this subsection is to 
study structures that do not deviate far from the uniform 
Salt and Pepper configuration discussed in the previous 
subsection. Because we have solved the uniform configu- 
ration exactly, the configurations which are not far from 
it are also treatable by the perturbation theory analysis 
i. e. expansion of the wire length (B) in terms of the 
deviation of densities of right and left eye neurons from 
the constant. This treatment determines which of the in- 
homogeneous phases (Stripe or Patch) is optimum. Also, 
comparison with the numerical results shows that the 
perturbation theory results hold even for big differences 
in density. 

We therefore consider a small repositioning of neurons, 
leading to the deviation of densities from constant Sn (r) . 
Because ur + nj, = uq 



n R (r) = n R + Sn (r) , 
n L (r) =fiL-8n (r) . 



(8) 



As this is only rearrangement the average of Sn (r) over 
the entire volume Sn (r) is zero, i. e. the total number 
of left and right eye neurons is not changed by the per- 
turbation. We then substitute these functions into our 
model (||) - (||) and calculate expansion of the wire length 
in the Taylor series in Sn (r) . It has the form: 

C = C SP + £d) + £(2) + . . . , 



(9) 



where £ SP is given by Eq. (Q), C^ oc Sn, C^ ex Sn 2 
are the first and the second order corrections to the wire 



length. From the condition Sn (r) = it follows that 
D- x > =0. The second order correction to the wire length 

is 



,SP 

Ik 



C^= I dvdv 1 J2 R l 

J i,k=R,L 

x {U lik (r - r') Sn, (r) Sn k (r') 
+ U 2 ik (r - r') Sn k (r) 5n k (r')} . 



(10) 



Here Sni (r) is the perturbation of density of neurons of 
i-th dominance (SriR = Sn, Snj, = —Sn,) and 



U U k(r)=0(R^-\r\) 



R SP 

n ik 



U 2ik (r) 



Mflgr 



Jdr"e(Rl k p -\r-r"\)9(R!£-\r"\) 



(11) 



(12) 



where 9{x) = 1, if x > 0, and 9(x) = 0, if x < 0. Because 
C^2ifc has the geometrical interpretation of the overlap be- 
tween two disks: 



U 2ik (r) 



1 

2^ 



2i& P 



4^ fc P 



1- 



2^ 



(13) 



2i& P - |r|) 



Using Eq. (||) we express the second order correction to 
the wire length (HG) as a pairwise density-density inter- 
action 

C {2) = J drdr'Sn (r) U (r - r') £n (r') , (14) 

where the "interaction potential" hi (r) is given by 

W(r)=Wi(r)+W 2 (r), (15) 



^1 = C/l_R_R + C/lLL — U\RL — U\LR, 



Ul — U2RR + U2LL + U2RL + U2LR- 



(16) 



(17) 



We notice that expression nj is similar to the Hamil- 
tonian used by Cowan and Friedman, 1991, which corre- 
sponds to Swindale, 1980, learning rules. The advantage 
of our approach is that we derive this expression from 
a single principle without assuming a particular form of 
"interaction potential" . In addition, we go on to solve 
this expression analytically. This allows us to map out a 
phase diagram which relates the ocular dominance pat- 
tern to biologically measurable connection rules without 
appealing to "Mexican hat" connection weights. 

We convert Eq.(|14|) into Fourier space using the prop- 
erty of a convolution 



£ (2) 



%U( q )\5n(cO\ 2 . 



(27T) 



(18) 



where hi (q) and Sh (q) are Fourier transforms of the "in- 
teraction potential" and the perturbation of density re- 
spectively. The Fourier transform of a function /(r) is 
defined as /(q) = f drf(r) exp(— iqr), where i = ^/~^l is 
the imaginary unity. Eq. (|18| ) is the central result of this 
subsection. 

Function hi (q) determines the changes in the total wire 
length due to the deviation of the neuron density from 
constant. For example, if the perturbation of density has 
the form of plane wave (5n = acos(qor)), the change in 
the total wire length is proportional to hi (qo) a 2 . Thus 
if hi (q) is negative at certain qo , such a perturbation de- 
creases the total wire length. It is therefore advantageous 
from the point of view of wire length economy to create a 



perturbation of density at this wave vector. In this case 
the uniform Salt and Pepper configuration (Sn = 0) is un- 
stable with respect to the formation of the OD patterns. 
Hence, negative function Li (q) indicates the formation 
an OD pattern. 

We therefore analyze the conditions at which the func- 
tion has negative values. Two statements can be made 
in this respect. First, assume that N s — N . Then 
hi (q) is never negative. This implies that Salt and Pep- 
per is optimum if N s — N . Second, consider hi (q) at 
Jr = II = 1/2 and arbitrary N s ^ N . In this case 
hi (q) always has negative values. This means that on the 
line of equal right-left eye occupancy Jr = Jl = 1/2 the 
OD patterns are always optimum, except for the point 
N s = N . We do not give the proofs of these properties 
due to the space limitations. 




FIG. 6. Function hi (q) calculated numerically for JV S = 10, 
No — 7, /b = 1/2, and no — 1. The value of wave vector 
corresponding to the most negative value of the function is 
denoted Qo- 



To illustrate these properties we show an example of 
hi (q) in Fig. g. The function obviously has negative 
values, signifying instability and an OD pattern forma- 
tion. The instability is strongest at the wave vector cor- 
responding to the most negative value of hi (q) . Indeed, 
creating the structure at this wave vector reduces the to- 
tal wire length most effectively. We predict therefore the 
spatial period of the OD pattern. For the case N s sa N 
shown in Fig. ^ (N s = 10 and iV = 9) function hi (q) 
reaches the most negative value at 



Qo 



Rrr 



ff SP 
n RL 



The spatial period of the OD pattern is therefore 
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(19) 



(20) 



In other words it is approximately equal to the diameter 
of the disc of connections. 

Competition between the Stripe and Patch phases Next 
we use the perturbation theory to calculate approxi- 
mately the wire lengths of different OD structures. Be- 
cause the structures are periodic the integral in Eq. JlS| ) 
can be reduced to the sum over the reciprocal lattice vec- 
tors Q: 



£(2) = iE^(Q)i^(Q)i 2 > 



(21) 



Q^O 



where A is the total area of the system. Different OD 
structures have different sets of Q and Sh(Q). For 
example, for Stripes Q x — 2im/ A, Q y = 0, where 
n = ±1, ±2, . . . and A is the spatial period of the struc- 
ture. The Fourier transform of density 
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For the triangular lattice of Patches Q x = Qq(1^/3/2), 
Q y = Qo(k + l/2), with l,k = ±1,±2, ... and Q = 
47r/Av / 3, where A is the lattice spacing. The Fourier 
transform of density 
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Based on Eqs. (pl|)-(|23]) we compare different OD struc- 
tures and generate the phase diagram similar to one given 
in the introduction (see Fig. [7]). 

0.6l 




FIG. 7. Ocular dominance phase diagram calculated in per- 
turbation theory. Range of the Stripe phase optimality is 
shown in black, Patch phase - in white, Salt and Pepper - in 
grey. 

Figures g and [7] have many similar features. First, the 
diagram is symmetric with respect to the line fn = 1/2. 
This is a consequence of the left-right eye symmetry of 
the general wire length functional (||) inherited by the 



second order functional (|1J). The reason for the exis- 
tence of such a symmetry is interchangeability of left and 
right eyes inherent to this model. If, for instance, in a 
given configuration one relabels left-eye neurons into the 
right-eye ones and wise versa, the wire length does not 
change. 

Second, the Salt and Pepper phase occupies a stripe 
around the line N s — N . The width of this line is 
given by \N S - N \ < 0.01N S . This is the result of the 
above-mentioned stability of Salt and Pepper on the line 
N s = N . As it is shown by the diagram the stability 
extends into some region around this line. 

Third, there is a transition between Stripes and 
Patches at fn fa 0.4 and 0.6. The region on the diagram 
corresponding to 0.4 < fn < 0.6 is almost completely 
occupied by the Stripes while the rest of the diagram 
{fn < 0.4 and fn > 0.6) by the Patches. We explain 
this in the framework of the perturbation theory. The 
main contribution to Eq. (|2lJ) comes from the terms with 
the smallest |Q|. This happens because both U (Q) and 
h (Q) decay very fast with the increase of |Q|. Stripes 
and Patches can approximately be compared using only 
the terms with the smallest |Q| = Qo- The two solutions 
have equal wire length if 



2U (Qo) |£«Stripes (Qo)| = 6W (Q ) |<5npatches (Qo)T 



(24) 



where the factors 2 and 6 are the numbers of the small- 
est wave length harmonics in the Stripe and Patch phases 
respectively. Using Eqs. ( f22[ ) and (g3|) to solve the lat- 
ter equation for fn we obtain numerically for the filling 
factor of the transition fn as 0.4. Due to the mentioned 
left-right eye symmetry of the model similar transition 
occurs at f R = l — 0.4 = 0.6. 

We would like to notice finally that comparison of the 
perturbation theory to exact calculations shows that the 
former works well even if the deviation of the density 
from constant is not small (~ 0.5no). Such a comparison 
shows that ([£-£ SP ] -£( 2 ))/£( 2 ) < 5%. In addition 
the perturbation theory provides a framework to under- 
stand numerous qualitative features of the phase diagram 
discussed above. 

Von der Malsburg, 1979, has surmised that there is a 
phase transition between Patches and Stripes driven by 
the cost of the left /right eye boundary. However he did 
not address different numbers of connections with same 
vs. other-eye neurons and made several different assump- 
tions (e.g. fixing the periodicity of the pattern). Thus 
our results offer a more complete description of the OD 
patterns while relying only on one principle - wire length 
minimization. 



III. DISCUSSION 

A. Comparison with experiment 

This theory relates functional requirements on a neu- 
ral circuit to its structural properties. In particular, the 
phase diagram relates the relative fractions of neurons, 
Jl, and of connections, N s /N , to the appearance of the 
OD pattern. Ideally, this theory could be tested by mea- 
suring these numbers experimentally and comparing the 
observed OD pattern to the one predicted by the phase 
diagram. However, we could not find data on the ratio 
N s /N and can only surmise that it is greater than one. 

We can partially test the theory by using the predic- 
tions of the phase diagram which are independent of the 
ratio N s /N . Fig.|l] shows that the transition from the 
Stripe and the L-Patch (R-Patch) phase takes place when 
h ~ O- 4 Ur ~ 0.4) for a wide range of N s /N . This 
number can be compared with the experimentally derived 
value of /t which is found from the relative area occupied 
by left-eye dominated neurons. The prediction that the 
Patch phase becomes optimal when one eye dominates 
is, indeed, non-trivial because there may be a system of 
alternating wide and narrow monocular stripes instead. 

We test this prediction on the data from macaque 
and Cebus monkey. The relative area occupied by the 
left /right eye depends on the location in VI. In para- 
foveal part of VI both eyes are represented equally, 
f R ss 0.5. In agreement with the phase diagram, the 
OD pattern consists of stripes. Farther from the fovea 
contralateral eye becomes dominant. The OD pattern 
becomes patchy there, just as expected from the phase 
diagram. We verify the location of the transition by using 
the following algorithm. We find /& for each point of the 
pattern by calculating the relative area occupied by the 
left/right regions in a window centered on that point and 
a few OD periods wide (dashed lines in Fig.||). Then we 
draw a contour corresponding to Jl = 0.4, FigJ|. Next 
we check visually whether the location of this contour is 
close to the transition from Stripes to Patches. Indeed, 
the large black contour in Fig.g coincides with the tran- 
sition indicating good agreement. 

In Cebus monkey the OD pattern has a similar transi- 



tion (Rosa et al., 1992). For monkey C06L from Rosa et 
al., 1992, we determine visually that along the horizon- 
tal meridian the transition occurs at the eccentricity of 
20 — 40 deg. According to the plot of the relative repre- 
sentations given in Rosa et al., 1992, /j, changes in the 
range 0.32 — 0.42 at these eccentricities. Our prediction 
of Jji — 0.4 falls into this interval. For the upper 45 de- 
gree meridian of the same monkey the transition occurs 
at the eccentricity of 30 — 40 degrees or at filling frac- 
tions 0.33 — 0.43. Again, the predicted value belongs to 
this interval. We conclude that this data agrees with our 
predictions although a more precise measurement would 



be helpful. 

In cats the OD patterns resemble Patches. In this case 
our theory implies that one eye should dominate. In fact, 
Shatz and Stryker, 1978, reported that the filling fraction 
of the contralateral eye in cat VI is greater than 0.5. This 
may explain the existence of Patches in cat VI. However, 



other authors (Anderson et al., 1988) claimed that both 
eyes are represented almost equally. More precise mea- 
surements of the ocular dominance are needed to make a 
conclusive judgment. 




FIG. 8. Transition between the Stripe and Patch phases oc- 
curs at theoretically predicted value /l. Shown is a fragment 
of the macaque ocular dominance pattern from Horton and 
Hocking, 1996. Neurons dominated by the left eye are grey 
and neurons dominated by the right eye are white. Black 
contours correspond to the value /l = 0.4 averaged over a 
window equal to the one shown (dashed lines). Transition 
from Stripe to Patch phase visually coincides with the black 
contour. 

B. Further Development of the Theory 

Next, we elaborate on several simplifying assumptions 
made in the paper. Although these assumptions should 
not affect our conclusions significantly, they are worth 
further exploration. 

First, the transition between Stripes and Patches may 
be more complex than discussed. We considered only 
two candidate phases: Stripes and a triangular lattice 
of circular Patches. It is possible that some intermediate 
phases become optimal near the transition. For example, 
Fig.|| hints that parallel chains of elongated Patches may 
give more efficient wiring. This would slightly modify our 
phase diagram. 

Second, we based our theory on looking for an opti- 
mal layout of neurons which minimizes total wire length. 



The considered structures are, therefore, regular and pe- 
riodic. However, developmental noise may lead to fluctu- 
ations in the OD pattern which reduce slightly its wiring 
efficiency. Although actual OD patterns contain such 
fluctuations we do not know whether these are due to 
suboptimal wiring or variations in the wiring rules from 
point to point. 

Different phases may have different stability in respect 
to noise. Judging from the data, the Stripe phase holds 
up well on the scale of a few periods. The Patch phase, 
however, does not show a regular triangular lattice. We 
think that this is because of a relatively small difference 
in wire length between the triangular and the square lat- 
tice of Patches. (It is about 0.5% of the total wire length, 
compared to 2% difference between Stripes and Patches 
for the upper left part of the phase diagram.) 

Third, we ignored the presence of numerous other cor- 
tical maps, such as orientation selectivity and retinotopy 
of receptive fields. This follows from the assumption that 
the mutual interaction (or coupling) between different 
maps is weak. For example, as mentioned in the Intro- 
duction, our motivation to neglect retinotopy comes from 
the magnitude of the receptive field scatter exceeding the 
width of ocular dominance stripes. The best justification 
for ignoring the coupling between the maps comes from 
the robustness of the observed OD patterns and the good 
agreement of our theory with experiment. 

Our theory can be expanded to address the interac- 
tion between different maps. Variables of other maps 
can enter the expression for the total wire length, EqJ|, 
through additional values of indices i and k, which so far 
reflect ocular dominance. Moreover, these indices can 
become continuous variables if the sums (g) are replaced 
by integrals. This would be appropriate for including in- 
teractions with retinotopic and orientational selectivity 
maps. 

Fourth, we applied our theory to the OD patterns as 
the best studied structure. Since our model is based on 
minimal assumptions, it can be applied to other binary 
structures such as cytochrome oxidase blobs. 

In conclusion, we explained the OD patterns in mam- 
malian VI by minimizing wire length given general func- 
tional considerations. Good agreement with experiment 
lends strong support to the notion that OD structures 
are adaptations to reduce intra-cortical wiring. 
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